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Abstract. The number of nucleon (NN) scaling of the directed flow vı and elliptic flow (v2) as well as the 
nuclear modification factor (Rep) are tested for light nuclei which are produced in 0.44 GeV 197 Au +° Au 
collisions at different impact parameters with two different in-medium nucleon-nucleon cross sections in a 
framework of an isospin-dependent quantum molecular dynamics (IQMD) model. In that energy domain, 
the emission of light nuclei can be well described by the squeeze-out phenonomenon. The results show a 
nice NN scaling behavior for flow parameters v1, v2 and Rep. These results demonstrate that the nucleon 
coalescence mechanism is responsible for nucleon-number scaling of above physical observables in squeeze- 
out region in heavy-ion collisions at intermediate energy. 
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1 Introduction 


In intermediate-energy heavy-ion collisions (HIC), the col- 
lective flow is one of the important observables which can 
reflect the early stage of heavy ion reaction, evolution dy- 
namics, in-medium nucleon-nucleon cross section [1-7], as 
well as the information on the nuclear equation of state 
(EOS) [8-16]. Many studies of anisotropic flows, namely 
the directed flow (vı) and the elliptic flow (v2) revealed 
the properties and origin of the collective motion through 
exploring the dependences for anisotropic flows on beam 
energy, fragment mass, impact parameter and so on in 
HIC [8-16]. 

The number-of-nucleon (NN) scaling of anisotropic flows 
was first proposed by Yan and Ma et al. in low-intermediate 
energy heavy ion collisions through a framework of a quan- 
tum molecular dynamics model [17,18]. This scaling states 
that when the anisotropic flows were scaled by the nucleon- 
number of light nuclei, these scaled flows will follow the 
same curve as a function of rapidity or transverse momen- 
tum per nucleon regardless of the mass of light nuclei [18]. 
This phenomenon can be naturally interpreted by a nu- 
cleonic coalescence mechanism for the formation of light 
nuclei at kinetic freeze-out stage in the reaction system. 
However, the earlier focus of the coalescence mechanism is 
only on the spectra of kinetic energy or momentum of light 
nuclei in heavy ion collisions [19] but not on flow behavior. 


Send offprint requests to: 
a Email: ygmaQsinap.ac.cn 


After the first prediction on number of nucleon scaling of 
flows in heavy-ion collisions at lower energies [17,18], the 
same concept of number-of-nucleon scaling of elliptic flow 
was followed by a dynamical model as well as a multiphase 
transport (AMPT) model calculations at ultra-relativistic 
energy [20-22]. Later on the number-of-nucleon scaling 
of anisotropic flows for light nuclei have been experimen- 
tally confirmed in Au + Au collisions at much higher en- 
ergies, such as 200, 62.4, 39, 27, 19.6, 11.5, and 7.7 GeV 
at the BNL Relativistic Heavy Ion Collider (RHIC) [23], 
which illustrates that the light nuclei at RHIC energies 
are also formed by nucleonic coalescence mechanism at 
the freeze-out stage but not directly by the quark coales- 
cence mechanism. At both lower and relativistic energies, 
the elliptic flow demonstrates the positive values: the for- 
mer stems from the collective rotational motion due to 
the dominated mean field, and the later origins from the 
strong pressure of the overlapping participant region. In 
both cases, the NN scaling has been examined. However, 
in middle energy region where the negative elliptic flow 
might emerge due to the squeeze-out behavior of particle 
emission, NN scaling of flows is not yet checked so far. 


On the other hand, the number-of-constituent-quark 
scaling of nuclear modification factor (NMF) for pions and 
protons as well as the number-of-nucleon scaling of NMF 
for light nuclei have been first proposed by one of the 
authors [24], and they were experimentally supported for 
Au + Au and Pb + Pb collisions at NA49, RHIC and 
LHC energies. However, the similar NN scaling behavior of 
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NMF is not yet examined in intermediate energy heavy ion 
collisions. In this energy domain, the isospin-dependent 
quantum molecular dynamics (IQMD) model [25-27] is 
very successful for studying reaction dynamics, such as 
for spectra analysis, flow analysis as well as the nuclear 
modification factor analysis [28]. Therefore it will be of 
very interesting to test if the Rep of light nuclei can be 
scaled by the number of nucleon scaling. 

Based upon the above arguments, in this article, we use 
an isospin-dependent molecular dynamics (IQMD) model 
to simulate 0.44 GeV Au + Au collisions at different 
impact parameters with two sets of in-medium nucleon- 
nucleon cross section. The anisotropic flows and Rep of 
different light nuclei especially in terms of the coalescence 


mechanism are systematically investigated, and the squeeze- 


out behavior was cleanly observed by the negative elliptic 
flow. In this work the Rep is obtained by dividing the spec- 
tra in the centrality of 0-10%, 10-20%, and 20-40% to the 
one in the centrality of 60-80%. Dependences of the mass 
number of the anisotropic flows and the nuclear modifica- 
tion factor are surveyed, and the nucleon-number scaling 
of light nuclei on the anisotropic flows and Rep are also 
found. 


2 Theoretical descriptions 
2.1 Definition of the directed and elliptic flows 


Anisotropic flow is defined as the different n-th harmonic 
coefficient v, of an Fourier expansion of the particle in- 
variant azimuthal distribution 


dN = 
rT «x1+2 5 Un cos(n¢), 
n=1 


(1) 
where the reaction plane was defined by a plane of the 
z-axis along the beam axis and x-axis along the impact 
parameter axis in the coordinate system and ¢ is the az- 
imuthal angle between the transverse momentum of the 
particle and the reaction plane. For a given rapidity the 
anisotropic flows are evaluated according to 


Un (pe) = (cos(ng)), (2) 


where (cos(n@)) denotes the average over the azimuthal 
distribution of particles with transverse momentum pz. 
Specifically, the directed flow 


and elliptic flow 


U2 = (cos(2¢)) = | ; 
Pt 


where py and py are components of momentum, and p; = 


\/p2 + pz is transverse momentum. 


2.2 Definition of nucleon-number scaling of nuclear 
modification factor 


The nuclear modification factor Rep is a probe to study 
the nuclear medium effect in nucleus-nucleus collisions, it 
is defined as a ratio of the particle yield in central collisions 
to that in peripheral collisions [29,24,28], i.e. 

central 
)| 


[d2 N/2rprdprdy/ (Nyin 
[2 N/2rprdprdy/ (Noin) P” 


(5) 


cp — 


where d?.N/2rprdprdy is the particle invariant differen- 
tial yield of the transverse momentum spectra pr, (Npin) 
means the average nucleon-nucleon binary collision num- 
ber per event. Now, we focus on the coalescence mecha- 
nism for the formation of light nuclei similar to that for 
the production [30] or elliptic flow [3,31] of hadrons in rel- 
ativistic heavy-ion collisions. The formation of light nuclei 
can be described by the coalescence mechanism as follows: 


(6) 


where P, = P, = Pa/A and A is the number of protons 
in light nuclei. The coefficient B4 is the probability for 
coalescence of nuclei with mass number A, which depends 
on momentum as well as the fireball volume in coordi- 
nate space [32]. The coefficient B4 can be extracted from 
data or calculated by the coalescence mechanism [32]. The 
number of nucleon scaling for the Rep of light nuclei can 
be deduced from Eq. 5 and Eq. 6 [24] : 


Re, = a (Rep (A+ pr))'/4 (P) 


Here, the index c and p stands for the central collisions and 
peripheral collisions, respectively. R¢, is the scaled nuclear 
modification factor to scale proton Rep together. Then we 
can try to scale Rep of light nuclei by the following formula 
[24]: 


c\ 1/A-1 


Ro (pr) = (Rep (A- pr))/* (es, 


and 


Bae ) as (9) 


it, (or) = Ry (wr) / (ZE 


2.3 The isospin-dependent quantum molecular 
dynamics model 


In the following discussion, we introduce the isospin-dependent 


quantum molecular dynamics (IQMD) model briefly. The 
quantum molecular dynamics model is a n-body transport 
theory, it describes heavy-ion reactions from intermediate 
to relativistic energies. The main parts of QMD trans- 
port model include the following: the initialization of the 
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projectile and the target, nucleon propagation in the ef- 
fective potential, the collisions between the nucleons in 
a nuclear medium, the Pauli blocking effect, and the nu- 
merical test. The isospin-dependent quantum molecular 
dynamics model is based on the QMD transport model 
with the isospin factors taken into account [25,33]. As we 
know, the main components of the dynamics in heavy- 
ion collisions (HICs) at intermediate energies include the 
mean field, two-body collisions, and Pauli blocking. There- 
fore, it is important for these three components to include 
isospin degrees of freedom in the IQMD transport model. 

In particular, the density-dependent Skyrme potential 
Usky reads when the momentum dependent potential is 
included 


Usry = a( L) + BCL)? + tain? e( £P + 14, 


£ (10) 
Po Po Po Po 


where p and p, are total nucleon density and its normal 
value, respectively. The parameters a, 6, y, t4 and € are 
related to the nuclear equation of state [34,35] and listed 
in Table I, where K = 200 or 380 MeV means the soft- or 
the stiff- momentum dependent potential, respectively. 


Table 1. The parameters of the interaction potentials. 


io B Y ta E K 
(MeV) | (MeV) (MeV) | (MeV) | (MeV) 
-390.1 320.3 1.14 1.57 21.54 200 
-129.2 59.4 2.09 1.57 21.54 380 


In this model, the particles are identified using a mod- 
ified isospin-independent coalescence description, i.e., two 
nucleons are assumed to share the same cluster if their 
centers are closer than a distance of 3.5 fm and their rel- 
ative momentum smaller than 0.3 GeV/c. If the nucleon 
is not bounded by any clusters, it is treated by an emit- 
ted (free) nucleon. In the IQMD model, the in-medium 
nucleon-nucleon cross section (NNCS) is represented by 
the formula: 


(11) 


where 7 is the in-medium NNCS reduction factor and 
a is the available experimental NNCS from Particle 
Data Group. The above in-medium NNCS reduction fac- 
tor ņ varied from 0 to 1 and the relationship between 7 
and the in-medium effect is opposite in this expression. 
In particular, the factor 7 ~ 0.2 has been found better to 
reproduce the flow data. 


3 Results and discussion 
Reactions for the collision system of 197 Au+1!97 Au at 0.4A 


GeV with different impact parameters have been simu- 
lated by IQMD model under the soft EOS with momentum- 


results are extracted at 200 fm/c when the reaction sys- 
tem has been at the freeze-out stage. Fig. 1 and Fig. 2 
shows the anisotropic flows vı as a function of rapidity 
and the flows per nucleon v;/A versus rapidity for test- 
ing the effect of the number of nucleon scaling, respec- 
tively, at different centralities with two sets of in-medium 
nucleon-nucleon cross section. In Fig. 1, we show the trend 
of the directed flow vı. We found that the directed flow 
(slope in mid-rapidity) is positive, indicating the repul- 
sive nucleon-nucleon interaction are dominant, which is 
opposite to negative slope vı at low energy [18] where 
the attractive mean field plays a key role. It also shows 
the absolute values of directed flow are larger for heavier 
clusters at a given rapidity. On the other hand, with the 
increase of NNCS reduction factor 7 from the upper pan- 
els to lower panels in Fig. 1, the directed flow values seem 
to drop, indicating the origin of flow most comes from 
the nucleon-nucleon collisions. After we perform the NN 
scaling of vı, we display Fig. 2 where the curves for dif- 
ferent clusters almost stay together around mid-rapidity 
by the nucleon-number scaled vı /A, which illustrates that 
the directed flow of the light clusters satisfies the number- 
of-nucleon scaling. 


The panels in Fig. 3 show transverse momentum de- 
pendence of elliptic flow v2 for light clusters at different 
centralities with two sets of in-medium nucleon-nucleon 
cross section. From the figure, it shows that elliptic flows 
are negative and decrease with the transverse momentum 
pt. The negative value of v2 reflects that the light clus- 
ters are preferentially emitted out of the reaction plane, 
i.e. so-called squeeze-out phenomenon, and particles with 
higher transverse momentum tend to be stronger nega- 
tive və values, which is opposite to the collision of lower 
beam energy [18]. In addition, the heavier the clusters, the 
larger the value of və. This behavior is similar to the di- 
rected flow. In intermediate energy, it is well known that 
both the attractive mean field and the repulsive nucleon- 
nucleon collisions play important roles. At lower energy, 
the mean field becomes dominant to contribute to the for- 
mation of a rotating compound system, so the positive el- 
liptic v2 is essentially induced by the collective rotational 
motion [13,36,37]. In that case, the elliptic flow stems from 
the attractive mean field. However, at mediate energy as 
here 0.4A GeV, the two-body collision gradually plays im- 
portant role and the collective expansion is therefore de- 
veloped. Although the transverse momentum dependence 
of və (the absolute values of v2) is similar to the previous 
results at RHIC energies, the mechanism of the flow is 
very different in the two different energy region. At RHIC 
energies, the elliptic flow is mainly driven by the stronger 
outward pressure [2], and here the elliptic flow is mainly 
due to the shadowing effects of overlapping participant 
zone. As for the NNCS effect, the larger ņ seems to reduce 
the absolute vg values as vı case. After we perform the 
similar scaling of nucleon number as v1, we observe a nice 
number of nucleon scaling exists for the elliptic flow per 
nucleon as a function of transverse momentum per nucleon 
as displayed in Fig. 4. In comparison to the number-of- 


dependent interaction (SM+MDI). In this work, the physics constituent quarks scaling of p;/n at RHIC energies, this 
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behavior looks apparently similar in Fig. 4. The results in- 
dicate that the production of light nuclei stems from the 
nucleonic coalescence mechanism. 


In Fig. 5 , we show the Rep of light nuclei (p,d,t) as 
a function of pr for different central-peripheral pairs in 
Au+Au collisions at E = 0.4A GeV. It is cleanly seen 
that the Rep curves for proton, deuteron, and triton are 
obvious different. The multiple nucleon-nucleon scattering 
effect [38] tends to transform the longitudinal momentum 
into the transverse momentum, and the effect becomes 
stronger with the increase in pr in HICs, which leads to 
larger Rep in the high pr region. On the other hand, in 
the low pr region, radial flow plays a major role in cen- 
tral collisions, which pushes protons to higher Rep regions 
and results in the smaller Rep at low pr. For protons, 
the strength of Rep enhancement is suppressed at high 
pr with the increase in impact parameter as well as in- 
medium NNCS factor. We noticed that this factor depen- 
dence of Rep is quite similar to the results of Rep obtained 
in Refs. [39,28] due to the multiple nucleon-nucleon scat- 
tering effect or/and the radial flow effect. 


The upper panels and lower panels in Fig. 5 shows the 
comparison of Rep of light nuclei with different 7 values. 
The 7 values at 0.0 and 0.5 were used in the simulation of 
Au + Au at 0.4A GeV collisions. In a previous work, the 
value 7 = 0.2, i.e., 80% of the free space nucleon-nucleon 
cross section was obtained [40]. In fact, the medium effect 
is different in various ranges of incident energy and mat- 
ter density [41]. The Rep of light nuclei has an increasing 
trend with pr which is explained by the multiple nucleon- 
nucleon scattering effect as well as radial flow [39], and its 
trend becomes more rapid for protons with pr in the low 
n value case because of the high collision rate between nu- 
cleons. Collisions become certainly less in higher 7 values 
and it makes the multiple nucleon-nucleon scattering ef- 
fect less significant. Based upon the above argument, Rep 
of protons is a good quantity for studying the effect of 
the in-medium NNCS, it means that we can investigate 
Rep to draw a conclusion on the in-medium effect by a 
quantitative comparison between the model and the data. 


After performance of the number of nucleon scaling for 
Rep, Fig. 6 shows the NN scaling of Rep for light nuclei, i.e., 


Re From this figure, through the number of nucleon scal- 
ing, it was found the values of Rep for deuteron and triton 
can be scaled to proton’s after using a constant factor (eg. 
for panels (a) and (d), a constant factor 1/0.76) on pro- 
ton’s Rep. These constant factors show no dependence on 
in-medium NNCS factor and weakly depend on centrality. 
All the above phenomena can be seen as the coalescence 
mechanism for formation of light nuclei. In comparison to 
the scaling behaviors of Rep for light nuclei produced from 
Pb + Pb at 17.2 GeV, Au + Au at 200 GeV, and Pb + Pb 
at 2.76 TeV [24], this behavior looks apparently similar 
in Fig. 6. The number-of-nucleon scaling of Rep is increas- 
ing as a function of pr, which could stem from the larger 
radial flow or multiple nucleon scattering for light nuclei 
as demonstrated in our earlier work [39]. Therefore it im- 
plies that the number of nucleon scaling of Rep for light 
nuclei supports the viewpoint of the coalescence mecha- 


nism for formation of light nuclei at the kinetic freeze-out 
stage, and this scaling phenomenon is the same as that 
for elliptic flow [17,18,23]. For the NNCS effect on Rip 
it shows insensitivity due to the cancel effect of the Rep 
ratio. 


4 Summary 


In summary, we have investigated the number of nucle- 
ons scaling of anisotropic flows, namely vı and v2 for light 
charged particles produced by +97 Au +197 Au at 0.4A GeV 
where the light particle emission has squeeze-out behavior 
in the framework of isospin-dependent quantum molecular 
dynamics model. The different centralities and in-medium 
NNCS factors are taken into account. The results show 
that the curves of both v,/A and v2/A for different light 
charged particles stay together, which means that there 
exists directed flow and elliptic flow scaling on the nucleon 
number for light charged particles. Then, the number of 
nucleons scaling of the nuclear modification factor Rep for 
light charged particles also have been investigated. The 
results demonstrate that the scaled Rep presents a nice 
overlap all together after considering a constant difference 
factor between light charged particles and proton’s, indi- 
cating a NN scaling of Rep is also satisfied. In light of 
present study, we propose experimental studies along this 
direction. 
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Fig. 1. (Color online) Rapidity dependence of vı for light nuclei for 0.44 GeV Au + Au collisions at different centralities. 
Squares represent for protons, circles for deuterons, up-triangles for tritons. From the left panel to right panel, it corresponds 


to different centralities of 0 - 10%, 10 - 20%, 20 - 40% and 40 - 60%. The upper row represents the calculations with a factor n 
of the in-medium nucleon-nucleon cross section 0, and the lower row with n = 0.5. 
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Fig. 2. (Color online) Same as Fig. 1 but for the vı/A as a function of rapidity. 
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Fig. 3. (Color online) The transverse momentum (p+) dependence of v2 for light nuclei in 0.4 A GeV Au + Au collision at 
different centralities (from left to right: 0 - 10%, 10 - 20%, 20 - 40%). Squares represent for protons, circles for deuterons, 
triangles for tritons. The upper row represents the calculations with a factor 7 of the in-medium nucleon-nucleon cross section 
0, and the lower row with 7 = 0.5. 
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Fig. 4. (Color online) Same as Fig. 3 but for v2/A versus p;/A. 


Fig. 5. (Color online) Rep for proton, deuteron and triton in Au + Au collisions at E = 0.4A GeV. From left to right, it 
corresponds to the ratio of different centralities of 0 - 10%, 10 - 20%, 20 - 40% to the centrality 40 - 60%. Two cases of the 
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in-medium NNCS are considered: 7 = 0.0 (upper row) and 0.5 (lower row). 
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Fig. 6. (Color online) Same as Fig. 5 but for the scaled Rep , i.e. Ree 


